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We carry out large-scale micromagnetic simulations which demonstrate that due to topological 
constraints, internal domain walls (Bloch lines) within extended domain walls are more robust than 
domain walls in nanowires. Thus, the possibility of spintronics applications based on their motion 
channeled along domain walls emerges. Internal domain walls are nucleated within domain walls 
in perpendicularly magnetized media concurrent with a Walker breakdown-like abrupt reduction of 
the domain wall velocity above a threshold driving force, and may also be generated within pinned, 
localized domain walls. We observe fast held and current driven internal domain wall dynamics 
without a Walker breakdown along pinned domain walls, originating from topological protection of 
the internal domain wall structure due to the surrounding out-of-plane domains. 

PACS numbers: 75.60.Ch, 75.78.Fg, 75.78.Cd 


During recent years, a lot of effort has been devoted to 
understand properties of magnetic domain walls (DWs) 
and their dynamics. A major driving force behind these 
studies is the emergence of next generation ICT compo¬ 
nents based on DWs, such as memory devices m and 
logic gates [3]. At the same time, magnetic DWs con¬ 
stitute a suitable playground to study several key funda¬ 
mental physics ideas and concepts, ranging from topology 
[4] to non-equilibrium critical phenomena [5j [6]. 

In general, DWs may have various internal degrees of 
freedom [7] which are essential for their magnetic held 
or spin-polarized current driven dynamics, and have re¬ 
cently been shown to be useful e.g. for channeling spin 
waves along extended DWs [S> In narrow ferromagnetic 
(nano)strips [both in soft Permalloy strips, and strips 
with a high perpendicular magnetic anisotropy (PMA)], 
the Walker breakdown [9], or the onset of precession 
of the DW internal magnetization m DW , leads to an 
abrupt decrease of the DW propagation velocity above 
the Walker held T>w or current density jw fT0HT2] . Such 
behaviour is captured also by simple one-dimensional 
(ID) models [12] of DWs in nanowires/strips. 

In wider PMA strips with longer DWs (such that the 
ID models are no longer applicable), one may expect that 
a Walker breakdown -like abrupt reduction of the DW 
propagation velocity still takes place, but that the related 
excitation(s) of m DW can no longer be spatially uniform 
D3I- Possible origins for the symmetry breaking leading 
to incoherent precession of m DW in different parts of the 
DW could be edge effects, and/or quenched disorder, in¬ 
teracting with the DW ill; [Ml 120] : these may include dis¬ 
locations, precipitates, grain boundaries, thickness fluc¬ 
tuations of the strip, etc. Here we explore the dynamics 
of extended DWs in wide CoPtCr PMA strips, with a 
Bloch wall equilibrium structure, using large-scale micro- 
magnetic simulations with and without quenched disor¬ 
der. We show that in wide enough strips and for driving 
forces exceeding a threshold value (T>w or jw)-> the in¬ 


ternal degrees of freedom of the DW are indeed excited 
inhomogeneously along the DW. These non-uniformities 
can be well-described as a set of internal DWs (IDWs), 
similar to Bloch lines observed e.g. in the context of 
bubble domains in garnet films [2TH25] . separating in¬ 
plane domains within the main DW, see Fig. [l] (d) for 
an example snapshot. Such IDWs, which resemble trans¬ 
verse DWs m found in narrow nanostrips with in-plane 
magnetization [27H29] , can also be nucleated within DWs 
pinned by notches, and subsequently driven along the 
localized main DW by magnetic fields or spin-polarized 



FIG. 1. (color online) Snapshots of field-driven DWs in per¬ 
fect CoPtCr strips with H ex t along z. (a) For a narrow strip 
(L y = 150 nm), the DW is a Bloch DW with m DW roughly 
along y for B ext < Bw, while (b) for H ex t > 2?w> m DW oscil¬ 
lates between spatially almost uniform configurations, (c) For 
a wider strip (L y = 1.2 /mi), m DW is uniform for H ex t < 5w, 
but (d) breaks into a dynamic pattern of in-plane domains 
separated by IDWs for B ext > 










2 


currents. We compute the driving force dependence of 
the velocity of the various IDW types observed, and find 
that they do not experience a Walker breakdown, due to 
topological constraints originating from the out-of-plane 
domains surrounding the IDWs. Thus, the possibility to 
use extended DWs as channels for fast IDW propagation 
emerges, with potential applications in spintronics. 

The micromagnetic simulations are performed using 
the GPU-accelerated micromagnetic code MuMax3 m v 
[32], solving the Landau-Lifshitz-Gilbert equation [33ll34h 

dm/dt = yHeff x m + am x dm/dt, ( 1 ) 

where m is the magnetization, 7 the gyromagnetic ra¬ 
tio, and H e ff the effective field, with contributions due 
to exchange, Zeeman, and demagnetizing energies. We 
consider CoPtCr strips of thickness L z = 20 nm [35] , 
and widths L y ranging from 150 nm to 3 /im, with the 
saturation magnetization M s = 3 • 10 5 A/m, exchange 
constant A = 10 -11 J/m, and the damping parameter 
a = 0.2. The first order uniaxial anisotropy constant 
K u = 2 • 10 5 J/m 3 , and an out-of-plane easy axis is con¬ 
sidered, to take into account the PMA nature of the strip 
|36 . The discretization cells have dimensions of = A y 
= 3 nm, and A z = 20 nm; we chose to use only one layer 
of computational cells in the ^-direction, since a finer dis¬ 
cretization in that direction did not change the results. 

We start by considering the field-driven dynamics of 
DWs in perfect CoPtCr strips, with two out-of-plane do¬ 
mains separated by a Bloch DW with the internal mag¬ 
netization m DW = -\-M s y as an initial state. We em¬ 
ploy a simulation window of length L x = 6 /im centered 
around and moving with the DW. When driven with a 
field of magnitude H ex t < B^(L y ) along z, we observe 
a slight, i? ex t-dependent tilting of the propagating DW, 
see Figs, [l] (a) and (c): as B e ^ t z rotates m DW away 
from the y direction, the DW tries to partially align it¬ 
self with its m DW , to minimize stray fields. Notice that 
this effect arises here without the Dzyaloshinskii-Moriya 
interaction m- For F> ext > F>w, the DW dynamics de¬ 
pends strongly on the strip width: for narrow strips [see 
Fig. 0(b)], m DW oscillates between spatially almost uni¬ 
form configurations, while for wider strips with longer 
DWs, m DW breaks into a dynamic, spatially non-uniform 
pattern, nucleated from the sample edges; snapshots of 
the magnetization configurations [Fig. |T] (d)] reveal that 
m DW is broken into a set of in-plane domains, separated 
by IDWs (Bloch lines). These patterns are highly dy¬ 
namic, with different parts of the DW either ahead or 
behind the average DW position. The IDWs move along 
the narrow (width ~ 9 nm) main DW, with annihila¬ 
tion and nucleation events of pairs of IDWs taking place 
repeatedly 1 38] : an example animation of this complex 
process is provided as Supplemental Material [39] . 

Fig. 1 (a) shows the steady state DW velocity udw as 
a function of H ex t for perfect strips of different widths 
L y . For all L^-values considered, udw exhibits the same 


FIG. 2. (color online) The average DW velocity vdw as a 
function of B ext in CoPtCr strips of different widths L y , in (a) 
perfect strips and (b) strips with a random 7 % standard de¬ 
viation variation of the anisotropy strength K u in each grain; 
an example of the grain structure is shown in the inset. Notice 
the strong size effect in the depinning field in (b). 


linear dependence on B ext for F> ex t < Byj(L y ), with 
B w (L y ) in the range 9-10 mT; above Fw, udw is sig- 
flcantly smaller than its (local) maximum obtained for 
smaller fields (indeed, a reduction of udw by a factor 
of a 2 is expected in the precence of many Bloch lines 
[231425] 5. and displays a relatively complex depedence on 
F> ex t, originating from the incoherent dynamics of m DW 
at different parts of the DW. Also current-driven DW dy¬ 
namics (not shown) exhibits a velocity drop and IDWs 
nucleated above jw m • 

In order to account for the effect of quenched disor¬ 
der, here assumed to originate from the poly crystalline 
structure of the strip [41], we construct grains using a 
Voronoi tessellation mm, with an average grain size 
of 11.9 nm [35], see the inset of Fig. [ 2 ] For each grain 
i, we assign a different, Gaussian distributed random 
anisotropy strength K Uj i, with mean K u and standard 
deviation a = 0.07 K u . Here, we consider strips of a fi¬ 
nite length L x = 6 /im. The DW is initialized at x = 1.5 
/im, with the origin being in the middle of the sample. 
The disorder results in a finite, cr-dependent depinning 
field, as evidenced by Fig. [ 2 ] (b), where averages over 
5 disorder realizations are presented. Notice the strong 
size effect in the depinning field, with the longer DWs 
depinning for smaller F> ex t- For longer DWs, m DW now 
precesses non-uniformly for all F> ex t values with udw > 0 , 
including those below Fw of the corresponding perfect 
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FIG. 3. (color online) IDWs may be generated within pinned DWs in disordered CoPtCr strips by field pulses. A field pulse of 
magnitude B ex t = 20 mT > Bw and duration of 5 ns, followed by relaxation in zero field, results in one pinned IDW, shown 
magnified in the last frame. The grey arrows indicate the direction of the in-plane DW magnetization around the IDW. 


system; thus, the non-uniform internal degrees of free¬ 
dom of the DW can also be induced by strong enough 
disorder, in addition to the sample edges. 

IDWs can also be created within localized, pinned 
main DWs; such setups could be useful to experimentally 
test our results. We consider here DWs pinned by trian¬ 
gular notches (with 50 nm long sides, and a 50 y/2 « 70.7 
nm long base) in strips of width L y = 1.2 /rm and length 
L x = 6 /im, with disorder a = 0.1 K u . By applying a 5 
ns long square field pulse of amplitude B ext = ±20 mT 
(i.e. above Bw) along z, the pinned DW bends as its 
central part propagates, and m DW exhibits spatially non- 
uniform dynamics (as B ext > Bw), see Fig. [3j After the 
pulse, the curved DW gets pinned by disorder, and may 
contain long-lived pinned non-uniformities in its m DW : 
stable 180° IDWs (Bloch lines) within the pinned main 
DW are formed with a roughly 40 % success rate (es¬ 
timated from an ensemble of 100 disorder realizations). 
Also more complex structures are sometimes observed, 
such as 360° IDWs. Fig. [3] shows a typical example 
of the process, with the end result of one head-to-head 
(H2H) 180° IDW, pinned by the disorder even after the 
relaxation time of 19 ns in zero field; the same process is 
also showed in a Supplemental Movie [39] . Magnification 
of the IDW magnetization (last frame of Fig. [3| reveals 
that the IDW resembles transverse DWs in narrow in¬ 
plane systems (e.g. Permalloy nanostrips) [26H28] . In 
addition to giving rise to different curvatures of the main 
DW, the direction (±z) of the applied field pulse affects 
the polarity of the IDWs created; B ext > 0 tends to lead 
to IDWs with > 0 (as e.g. in Fig. [3|, whereas 

B ex t < 0 gives mostly rise to IDWs with mJP w < 0; 
for both polarities, the IDW may have either a H2H or 
tail-to-tail (T2T) configuration. 

Examples of these four IDW configurations are pre¬ 
sented in Fig. [4j Due to the geometry shown in Fig. 
[4j the IDWs can be driven by applying fields along ±y 


directions. To study the velocity-field characteristics of 
the IDWs, we first move all of them from their random 
initial positions to close to the lower edge of the strip 
by applying field pulses of a small magnitude. Then, a 
driving field B ext is applied along +y or —y direction, 
to move the IDW towards the upper edge; the sign de¬ 
pends on whether the IDW has a H2H or T2T structure. 
The resulting GDw(^ext) curves are shown in Fig. [5j 
separately for the four possible IDW configurations, re¬ 
vealing that for B ext above a small depinning field, the 
T2T IDWs move significantly faster than the correspond¬ 
ing H2H IDWs. This can be understood by considering 
the force due to B ext on the IDW (with m IDW along ±x) 
in directions perpendicular to the main DW, shown with 
white arrows in Fig. [4j for H2H IDWs, the force direction 
coincides with that due to the line tension of the curved 
main DW [Fig. [4] (a) and (b)]. Thus, during the dynam¬ 
ics of the IDW, part of the energy of the driving field is 
dissipated in partial straightening of the main DW, re¬ 
sulting in a lower uidw for H2H IDWs. In contrast, for 
T2T IDWs, the perpendicular force on the IDW due to 
the driving field points in the opposite direction from that 
of the curvature-induced force: thus, less motion of the 
main DW takes place during the IDW dynamics, and a 
larger fraction of the field energy is available to move the 
T2T IDW, resulting in a larger uidw- This is confirmed 
by movies provided as Supplemental Material [39] : the 
motion of T2T IDWs is noticeably smoother than that of 
H2H IDWs. As a further check, we consider also a pure 
system with an artificially generated straight main DW 
containing an IDW; the resulting uiDw(^ext) curves for 
the four possible configurations are shown as solid black 
lines in Fig. |5j All cases exhibit the same niDw(^ext) 
behaviour, which is intermediate between those found 
for curved main DWs with H2H and T2T IDWs, respec¬ 
tively. Thus, the perpendicular forces still play a role in 
energy dissipation, but less than in the curved H2H case. 
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FIG. 4. Snapshots of IDWs of different structures, driven from 
bottom to top by a field B ex t along the d -y direction, shown 
as black arrows. The grayscale indicates the magnitude of 
m x . IDWs in (a) and (b) are of the H2H type, whereas (c) 
and (d) have T2T configurations. Depending on the sign of 
m IDW , F> ext gives rise to a force acting on the internal DW 
either towards left or right, as shown by the white arrows. 



B ext t mT ] 


FIG. 5. (color online) The average (over 4 disorder realiza¬ 
tions) IDW velocity i’idw as a function of F> ex t for the four 
possible IDW configurations. The black lines show the cor¬ 
responding data for perfect strips with artificially generated 
straight main DWs, each containing an IDW of one of the 
four different kinds (H2T,T2T, m™ w > 0, m™ w < 0). 

The ^iDw(^ext) curves in Fig. [5] exhibit two note¬ 
worthy features: (i) uidw grows sub-linearly with T? ext , 
indicating that an increasing fraction of the field energy 
is dissipated in other processes than IDW motion when 
T> ext is increased, (ii) No clear velocity drop, a signa¬ 
ture of Walker breakdown, is observable; this is the case 
even if we consider relatively large fields up to 15 mT. In¬ 
spection of the magnetization dynamics of the IDWs (see 
Supplemental Material [39]) reveals that the direction of 
m IDW is indeed preserved during the dynamics. This can 
be understood to follow from the peculiar topology of the 
system at hand: unlike nanostrips, the main DW along 
which the IDW is propagating does not have free bound¬ 
aries; these are crucial e.g. for nucleation of antivortices, 
mediating the magnetization reversal of transverse DWs 
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FIG. 6. (color online) The IDW velocity 'Uidw as a function of 
the current density j, for different values of the non-adibatic 
parameter £, in a narrow strip with the main DW along the 
long axis of the strip; part of such a strip is shown in the inset. 
For |j|-values larger than those shown, the domain structure 
breaks down. 


in in-plane nanostrips. Instead, here the IDWs are sur¬ 
rounded by the two out-of-plane domains, which topo¬ 
logically protect m IDW from being flipped, leading to the 
absence of Walker breakdown, and consequently to rela¬ 
tively large field-driven IDW velocities (up to 150 m/s for 
T2T IDWs). A similar absence of the Walker breakdown 
of topological origin has been reported in ferromagnetic 
nanotubes [42] . 

Finally, we consider current-driven IDW dynamics in 
narrow perfect CoPtCr strips with the main DW located 
in the middle of the strip and oriented along its long axis 
(inset of Fig. [6|. By extending Eq. 0 with spin-transfer 
torque terms [43] , we observe a simple linear dependence 
of uidw on the current density j for £ = a and £ = 2a 
(with £ the non-adiabatic parameter) up to uidw = 490 
m/s for £ = 2cq whereas for £ = 0, uidw = 0 for all j (Fig. 
|6|. Thus, Walker breakdown is absent also in current- 
driven IDW dynamics, due to the topological protection 
discussed above; this implies also that for £ = 0, IDWs 
remain intrinsically pinned m for any j. 

To conclude, we have shown that for DWs in wide 
PMA strips, a Walker breakdown -like abrupt reduction 
of the DW propagation velocity is concurrent with nucle¬ 
ation of internal in-plane domains separated by internal 
DWs (Bloch lines), resulting in a hierarchical DW struc¬ 
ture, with DWs within DWs. The IDWs can be driven 
along the main DW, and the absence of Walker break¬ 
down in their dynamics could lead to interesting possibil¬ 
ities for spintronics applications where DWs would serve 
as guides for fast IDW propagation. 

We thank Mikko J. Alava and Jan Vogel for useful com¬ 
ments. This work has been supported by the Academy 
of Finland through its Centres of Excellence Programme 
(2012-2017) under project no. 251748, and an Academy 
Research Fellowship (LL, project no. 268302). We ac¬ 
knowledge the computational resources provided by the 
Aalto University School of Science “Science-IT” project, 

























5 


as well as those provided by CSC (Finland). 


* touko.herranen@aalto.fi 

[1] S. S. Parkin, M. Hayashi, and L. Thomas, Science 320, 
190 (2008). 

[2] S. Parkin and S.-H. Yang, Nat. Nanotechnol. 10, 195 
(2015). 

[3] D. A. Allwood, G. Xiong, C. Faulkner, D. Atkinson, 
D. Petit, and R. Cowburn, Science 309, 1688 (2005). 

[4] O. Tchernyshyov and G.-W. Chern, Phys. Rev. Lett. 95, 
197204 (2005) 

[5] S. Zapperi, P. Cizeau, G. Durin, and H. E. Stanley, Phys. 
Rev. B 58, 6353 (1998). 

[6] S. Papanikolaou, F. Bohn, R. L. Sommer, G. Durin, 
S. Zapperi, and J. P. Sethna, Nature Phys. 7, 316 (2011). 

[7] V. Lecomte, S. E. Barnes, J.-P. Eckmann, and T. Gia- 
marchi, Phys. Rev. B 80, 054413 (2009) 

[8] F. Garcia-Sanchez, P. Borys, R. Soucaille, J.-P. Adam, 
R. L. Stamps, and J.-V. Kim, Phys. Rev. Lett. 114, 
247206 (2015). 

[9] N. L. Schryer and L. R. Walker, J. Appl. Phys. 45, 5406 
(1974). 

[10] E. Martinez, J. Phys. Condens. Matter 24, 024206 

( 2012 ). 

[11] E. Martinez, Adv. Cond. Matter. Phys. 2012 (2012). 

[12] A. Thiaville and Y. Nakatani, in Spin dynamics in con¬ 
fined magnetic structures III (Springer, 2006) pp. 161— 
205. 

[13] J. Zebrowski and A. Sukiennicki, J. Appl. Phys. 52, 4176 
(1981). 

[14] J. Leliaert, B. Van de Wiele, J. Vandermeulen, A. Coene, 
A. Vansteenkiste, L. Laurson, G. Durin, B. Van Waeyen- 
berge, and L. Dupre, Appl. Phys. Lett. 106, 202401 
(2015). 

[15] J. Leliaert, B. Van de Wiele, A. Vansteenkiste, L. Laur¬ 
son, G. Durin, L. Dupre, and B. Van Waeyenberge, J. 
Appl. Phys. 115, 233903 (2014). 

[16] B. Van de Wiele, L. Laurson, and G. Durin, Phys. Rev. 
B 86, 144415 (2012) 

[17] Y. Nakatani, A. Thiaville, and J. Miltat, Nature Mater. 
2, 521 (2003). 

[18] H. Min, R. D. McMichael, M. J. Donahue, J. Miltat, and 
M. Stiles, Phys. Rev. Lett. 104, 217201 (2010). 

[19] J. Leliaert, B. Van de Wiele, A. Vansteenkiste, L. Laur¬ 
son, G. Durin, L. Dupre, and B. Van Waeyenberge, Phys. 
Rev. B 89, 064419 (2014) 

[20] X. Jiang, L. Thomas, R. Moriya, M. Hayashi, 


B. Bergman, C. Rettner, and S. S. Parkin, Nat. Com- 
mun. 1, 25 (2010). 

[21] V. Volkov and V. Bokov, Phys. Solid State 50, 199 
(2008). 

[22] S. Konishi, IEEE Trans. Magn. 19, 1838 (1983). 

[23] A. Malozemoff and J. Slonczewski, Phys. Rev. Lett. 29, 
952 (1972). 

[24] J. Slonczewski, J. Appl. Phys. 45, 2705 (1974). 

[25] A. Malozemoff and J. Slonczewski, Magnetic domain 
walls in bubble materials (Academic Press, 1979). 

[26] R. Hertel and A. Kakay, J. Magn. Magn. Mater. 379, 45 
(2015). 

[27] R. D. McMichael and M. J. Donahue, IEEE Trans. Magn. 
33, 4167 (1997). 

[28] Y. Nakatani, A. Thiaville, and J. Miltat, J. Magn. Magn. 
Mater. 290, 750 (2005). 

[29] V. Estevez and L. Laurson, Phys. Rev. B 91, 054407 
(2015). 

[30] A. Vansteenkiste, http://mumax.github.io 

[31] A. Vansteenkiste and B. Van de Wiele, J. Magn. Magn. 
Mater. 323, 2585 (2011). 

[32] A. Vansteenkiste, J. Leliaert, M. Dvornik, M. Helsen, 
F. Garcia-Sanchez, and B. Van Waeyenberge, AIP Adv. 
4, 107133 (2014). 

[33] T. L. Gilbert, IEEE Trans. Magn. 40, 3443 (2004). 

[34] W. F. Brown, Micromagnetics , 18 (Interscience Publish¬ 
ers, 1963). 

[35] Y. Kubota, L. Folks, and E. E. Marinero, J. Appl. Phys. 
84, 6202 (1998). 

[36] D. Weller, A. Moser, L. Folks, M. E. Best, W. Lee, M. F. 
Toney, M. Schwickert, J.-U. Thiele, and M. F. Doerner, 
IEEE Trans. Magn. 36, 10 (2000). 

[37] O. Boulle, S. Rohart, L. Buda-Prejbeanu, E. Jue, 
I. Miron, S. Pizzini, J. Vogel, G. Gaudin, and A. Thiav¬ 
ille, Phys. Rev. Lett. Ill, 217203 (2013). 

[38] A. Zvezdin, A. Popkov, and I. Yarema, Zh. Eksp. Teor. 
Fiz 98, 1070 (1990). 

[39] See Supplemental Material at [URL will be inserted by 
publisher] for movies illustrating the details of the do¬ 
main wall dynamics discussed in the paper. 

[40] J. Iwasaki and N. Nagaosa, J. Phys. Soc. Jpn. 84, 083701 
(2015). 

[41] We expect other sources of disorder, leading to a random 
energy landscape interacting with the DW, to give rise 
to similar results. 

[42] M. Yan, C. Andreas, A. Kakay, F. Garcia-Sanchez, and 
R. Hertel, Appl. Phys. Lett. 99, 122505 (2011). 

[43] S. Zhang and Z. Li, Phys. Rev. Lett. 93, 127204 (2004). 

[44] T. Koyama, D. Chiba, K. Ueda, K. Kondou, H. Tani- 
gawa, S. Fukami, T. Suzuki, N. Ohshima, N. Ishiwata, 
Y. Nakatani, et a/., Nat. Mater. 10, 194 (2011). 



